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Un choix de dix thèmes étudiés cette année pour revoir les principes de l'algorithmique. 


Suites : limite infinie 

Suites : suite convergente 

Suites : somme des termes 

Équations : résolution par balayage 

Équations : résolution par dichotomie 

Calcul intégral : méthode des rectangles 

Calcul intégral : méthode de Monte-Carlo 

. Probabilités : simulation d'une expérience aléatoire 
Probabilités conditionnelles : évolution d'un processus 
. Lois continues : approximation par la loi normale 
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Algorithmique n° 1 Suites 


Limite infinie d’une suite 


On étudie la suite définie pour tout n > 0 par 


nm 
ue = 5 + (1) 


Il s’agit d'examiner son comportement lorsque n devient grand et plus précisément d'illustrer la définition 
suivante du cours : 


on dit que la suite (un) tend vers +o lorsque pour tout réel K il existe un rang n à partir duquel 
les termes de la suite sont supérieurs ou égaux à K 


Dans ce but on considère l’algorithme décrit par le pseudocode suivant ; 


VARIABLES 
n, u, p sont des nombres 
N est un nombre 

DÉBUT 
p prend la valeur 1 
N prend la valeur 9 
POUR n allant de O à N 


n 
u prend la valeur 2 + D 


afficher u 
p prend la valeur -1 * p 
FIN de la boucle POUR 
FIN 


Questions 


1. Faire fonctionner l’algorithme sur une feuille de papier de façon à obtenir les 10 premiers termes de 
la suite (u,) . Pour cela on peut reproduire et compléter le tableau suivant : 


n Ro A rez 
u = 
p 1 
N || 9 


2. Taper et tester l'algorithme avec Algobox : vérifier que les valeurs affichées sont les mêmes que celles 
qui sont dans le tableau de la question précédente. 
Quelle modification faut-il apporter pour afficher les 20 premiers termes de la suite ? 
3. On veut maintenant adapter l’algorithme de façon à afficher le premier entier n tel que un > 500 
a) Pourquoi une boucle POUR n’est pas adaptée à ce genre de question ? 
b) Modifier l'algorithme précédent en remplaçant la boucle POUR par une boucle TANT QUE 
- la condition devra être : u < 500 
- la variable N n’a plus aucune utilité 
c) Taper et tester l’algorithme modidié. 
4. Démontrer que pour tout n > 1002 on a bien u, > 500 
(le résultat est admis pour un réel À quelconque) 


Algorithmique n° 2 Suites 


Suite convergente 


On étudie la suite définie pour tout n > 0 par 


3 — Un 
5) 

Il s’agit d'examiner son comportement lorsque n devient grand et plus précisément d'illustrer la définition 

suivante du cours : 
on dit que la suite (un) tend vers le réel L lorsque tout intervalle ouvert contenant L contient tous les 
termes de la suite à partir d’un certain rang 


UQ — 28 et Un+1 = 


Dans ce but on considère l’algorithme décrit par le pseudocode suivant ; 


VARIABLES 
n, u sont des nombres 
N est un nombre 
DÉBUT 
u prend la valeur 28 
N prend la valeur ...... 
POUR n allant de 1 à N 
u prend la valeur ............... 
afficher u 
FIN de la boucle POUR 
FIN 


Questions 


1. Compléter le pseudocode de façon que l’algorithme affiche les 15 premiers termes de la suite. 
2. Taper et tester l’algorithme avec Algobox. 
Quelle conjecture peut-on faire à propos de la limite éventuelle de la suite ? 
3. On veut maintenant adapter l’algorithme de façon à afficher le premier entier n tel que u, soit dans 
l'intervalle ]0, 4 ; 0,6[ : on doit remplacer la boucle POUR par une boucle TANT QUE. 
a) La condition « u est dans l'intervalle 10,4 ;0,6[ » se traduit par : [0.4 < u ET u < 0,6 
Quel est son contraire logique ? 
b) Modifier l'algorithme précédent en utilisant une boucle TANT _ QUE 
c) Taper et tester l’algorithme modidié. 


4. Déterminer le premier rang no pour lequel un € ]0,49 ; 0,51 
Démontrer ensuite que si un € ]0,49 ;0,51{, alors un+1 € ]0,49 ;0,51[ 
Quel résultat obtient-on pour la suite (un) ? 

5. Le résultat précédent est admis dans le cas général : tout intervalle ouvert contenant 0,5 va contenir 
tous les termes de la suite à partir d’un certain rang. 


Algorithmique n° 3 Suites 


Somme des termes d’une suite 


Le Lièvre est à la poursuite de la Tortue. 
L T 


| l 
T T > 


La tortue a encore 5 m d'avance mais elle s’est endormie. 

Voyant cela, le lièvre ralentit son allure. La première seconde, il fait un bond de 1m, la seconde 
suivante il fait un bond de 0,50 m et ainsi de suite : chaque seconde, il fait un bond moitié moindre 
que le bond précédent. 

Au bout de combien de secondes le lièvre aura-t-il rattrapé la tortue ? 


Mise en forme du problème. 

On note u, la distance en mètres que le lièvre à parcouru au bout de n secondes. 

Puisque la première seconde le lièvre bondit de 1m , on aura : u1 = 1 

Sachant que la seconde suivante, il fait un bond de 0,50m, on aura : u2 = 1 +0,50 = 1,50 
Et en général, on pose | uy = 1 + 0,50 + 0,25 + etc... | (la somme comprenant n termes) 


L’algorithme suivant permet de calculer la distance parcourue au bout de n secondes. 


VARIABLES 
i et b sont des nombres 
n et u sont des nombres 
DÉBUT 
u prend la valeur 0 
b prend la valeur 1 
n prend la valeur ...... 
POUR i allant de 1 àn 
u prend la valeur ............... 
b prend la valeur ............... 
FIN de la boucle POUR 
afficher u 


FIN 
Questions 
1. Compléter le pseudocode de façon que l’algorithme affiche la distance parcourue au bout de 5 secondes 
en tenant compte des informations suivantes : 
+ n est le nombre de secondes au bout desquelles on calcule la distance parcourue 
x i est un compteur 
x u est la distance parcourue au bout de n secondes 
x b est le bond fait par le lièvre : il est donc divisé par 2 chaque seconde 
2. Taper et tester l’algorithme avec Algobox. 
3. Modifier l'algorithme de façon que l’utilisateur puisse rentrer le nombre de secondes au bout desquelles 
il veut connaître la distance parcourue par le lièvre. 
4. Proposer une réponse à la question initiale en testant plusieurs valeurs de n 
5. Démontrer la conjecture précédente 


(on pourra utiliser le résultat du cours sur la somme des termes d’une suite géométrique pour obtenir 
un majorant de u, et comparer ce nombre avec 5) 


Algorithmique n° 4 Fonctions 


Résolution approchée d’une équation par balayage 


Figure. 


On a représenté ci-dessus la courbe d’équation y = cos x et la droite d’équation y = x. 
On veut déterminer l’abscisse du point d’intersection des deux courbes qui est visible sur cette figure. 
On peut observer graphiquement que cette abscisse est comprise entre 0 et 1. 


Mise en forme du problème. 

On note f(x) la différence entre les deux fonctions représentées : 
f(x) = x — cosx 

Il est alors équivalent de résoudre l’équation | f(x) = 0 

Ne sachant pas déterminer une solution exacte, on se contente d’un encadrement de celle-ci. 


L’algorithme suivant teste si la fonction f change de signe entre deux abscisses x1 et x2 de l'intervalle 
[0 ; 1]. Une série de 10 tests est effectuée. 


VARIABLES 
a , b , p , xi , x2 , yl et y2 sont des nombres 
i est un nombre 
DÉBUT 
a prend la valeur 0 
b prend la valeur 1 
p prend la valeur (b - a) / 10 
xi prend la valeur a 
x2 prend la valeur xi + p 
POUR i allant de 1 à 10 
yi prend la valeur xi - cos(x1) 
y2 prend la valeur ............... 
si yi*xy2 < O afficher x1 et x2 
xi prend la valeur xi + p 
x2 prend la valeur ............... 
FIN de la boucle POUR 
FIN 


Questions 


1. Calculer f’(x) et démontrer que cette dérivée est strictement positive sur [0 ; 1] 
Démontrer que l'équation f(x) = 0 a une solution a et une seule sur cet intervalle. 
2. Compléter le pseudocode de façon que l’algorithme affiche les bornes de l'intervalle qui contient @ 
en tenant compte des informations suivantes : 
x p est le pas dont on augmente x1 et æ2 à chaque itération 
x i est un compteur 


ce 


Taper et tester l’algorithme avec Algobox. Indiquer l'encadrement de a qui est obtenu. 
4. On veut améliorer la précision de la réponse. 
Pour cela, on va exécuter trois fois la même boucle POUR comme le suggère le pseudocode suivant : 


VARIABLES 
a , b , p , xi , x2 , yl et y2 sont des nombres 
i et j sont des nombres 
DÉBUT 
a prend la valeur 0 
b prend la valeur 1 
POUR j allant de 1 à 3 
p prend la valeur ............... 
xi prend la valeur ............... 
x2 prend la valeur ............... 
POUR i allant de 1 à 10 
yi prend la valeur ............... 
y2 prend la valeur ............... 


SI yixy2 < O ALORS 
afficher xl et x2 
a prend la valeur xi 
b prend la valeur x2 
FIN de l’instruction SI 


xi prend la valeur xi + p 
x2 prend la valeur ............... 
FIN de la boucle POUR 
FIN de la boucle POUR 
FIN 


5. Compléter et tester l’algorithme proposé. 
6. Comment l’adapter de façon que l’utilisateur puisse donner au départ la précision souhaitée ? 


Algorithmique n° 5 Fonctions 


Résolution approchée d’une équation par dichotomie 


On sait que la fonction exponentielle est continue et strictement croissante et que ses limites en —c et 
en +oo sont respectivement 0 et +oo . Par conséquent l’équation 


er =? 


a une solution a et une seule qu’on se propose d’encadrer par « dichotomie » . 


Il est immédiat de vérifier que a est entre 0 et 1 : on a donc au départ un encadrement d’amplitude 1 . 
On va tester le milieu de cet intervalle, c’est-à-dire 0, 5 . Il y a deux possiblités : 

- si e0 > 2 alors a est entre 0 et 0,5 

- sie < 2 alors à est entre 0,5 et 1 
On obtient ainsi un encadrement d'amplitude 0,5. 


L’algorithme suivant applique cette méthode jusqu’ à obtenir une précision suffisante : 


VARIABLES 
X1 , x2 , xm , yi , y2 et ym sont des nombres 
p est un nombre 
DÉBUT 
xi prend la valeur 0 
x2 prend la valeur 1 
p prend la valeur x2 - xl 
TANT QUE p > 0.1 
xm prend la valeur (xi + x2) / 2 
ym prend la valeur exp(xm) 
si ym < 2 alors x1i prend la valeur xm 
SION: aan euece ur dort 
p prend la valeur ............... 
FIN de la boucle TANT QUE 


Questions 


Compléter l’algorithme de façon qu’il affiche les bornes d’un intervalle d'amplitude inférieure à 0,1. 
Taper et tester l’algorithme avec Algobox. Indiquer l'encadrement de a qui est obtenu. 

Comment améliorer la précision du résultat ? 

Comment l’adapter de façon que l’utilisateur puisse donner au départ la précision souhaitée ? 


One 2 EU 


On considère l'équation —2° — 5x + 12 = 0. 
Démontrer que cette équation a une solution et une seule dans R. 
Adapter l’algorithme précédent de façon à obtenir un encadrement de cette solution d'amplitude 


inférieure à 10° . 
Combien d’itérations de la boucle TANT _ QUE ont été nécessaires ? 


Algorithmique n° 6 


Calcul d’aire 


Approximation d’une aire par la méthode des rectangles 


On souhaite calculer l’aire sous la courbe de la fonction inverse entre 1 et 2. 


Figure. 


+ 
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L’intervalle [1 ; 2] à été partagé en 5 petits intervalles d'amplitude 0,2 . 
Le domaine dont on veut calculer l’aire est contenu dans la réunion de 5 petits rectangles dits « à 
gauche » de hauteurs respectives f(1) , f(1,2) , f(1,4), f(1,6) et f(1,8). 

Ce même domaine contient la réunion de 5 petits rectangles dits « à droite » de hauteurs respectives 


f(,2) , f(14), (1,6), f(1,8) et (2) . 


Questions 


1. Déterminer avec la calculatrice l’aire de la réunion des rectangles « à gauche » puis l’aire de la 
réunion des rectangles « à droite ». 


En déduire un encadrement de l’aire du domaine. Comparer avec la valeur exacte. 


2. On veut améliorer la précision du résultat. 
On note n le nombre d’intervalles construits entre 1 et 2 et h le pas utilisé . 
Chaque rectangle a pour base un intervalle {a ; b] dont les valeurs successives sont [1 ; 1,2], [1,2 ;1,4] 


éte…:.11,8:2). 


Dans un premier temps. on considère l'algorithme suivant 
? 


i ,n ,h ,a 


3 b 2 


G et D sont des nombres 


D D 55 


prend la valeur 5 
prend la valeur 1/n 
prend la valeur 1 
prend la valeur a +h 


POUR i allant 
G prend la 
D prend la 
a prend la 
b prend la 
FIN de POUR 
afficher G 
afficher D 


de 1 à n 


valeur 
valeur 


valeur a +h 
valeur a +h 


Compléter cet algorithme de telle façon qu'il affiche en sortie la somme des aires des rectangles à 
gauche et la somme des aires des rectangles à droite. 

3. On souhaite maintenant que l'utilisateur indique en entrée le nombre n d’intervalles à utiliser. 
Proposer et tester un algorithme qui réponde à cette question. 

4. Une assez bonne approximation de l’aire est obtenue en faisant la moyenne des deux approximations 
précédentes (méthode des trapèzes) : modifier et tester l’algorithme précédent de façon à implémenter 
cette méthode. 


1 
5. Déterminer une valeur approchée de | 4 eŸ dé. Quelle est la valeur exacte ? 
0 


Corrigé Calcul d’aire 


1. Aire des rectangles à gauche : G = [f(1) + f(1,2) + f(1,4) + f(1,6) + f(1,8)] x 0,2 & 0, 746 
Aire des rectangles à droite : D = [f(1,2) + f(1,4) + f(1,6) + f(1,8) + f(2)] x 0,2 & 0,646 
Encadrement de l’aire sous la courbe : 0,646 < À < 0,746 


2 
Valeur exacte : À = | _ = [In +]? — In2 & 0,693 
1 
2. Algorithme permettant de calculer l’aire des rectangles à gauche et l’aire des rectangles à droite. 
1 ,n ,h ,a , b , G et D sont des nombres 
n prend la valeur 5 
h prend la valeur 1/n 
a prend la valeur 1 
b prend la valeur a +h 
POUR i allant de 1 à n 
G prend la valeur G+h/a 
D prend la valeur D +h /Db 
a prend la valeur a +h 
b prend la valeur b +h 
FIN de POUR 
afficher G 
afficher D 
3. Algorithme avec choix du nombre d’intervalles. 
1 ,n ,h ,a , b , G et D sont des nombres 


etc ......... 
4. Implémentation de la méthode des trapèzes. 
1 ,n ,h,a ,b,T ,G et D sont des nombres 
lire n 
h prend la valeur 1/n 
a prend la valeur 1 
b prend la valeur a + h 
POUR i allant de 1 à n 


ELEC: as aautets 

FIN de POUR 

T prend la valeur (G + D) / 2 

Afficher T 

Résultats numériques : 

n 5! 10 510 
G 0.74563492 0.7187714 0.69817218 
D 0.64563492 0.6687714 0.68817218 
T 0.69563492 0.6937714 0.69317218 


Algorithmique n° 7 Calcul d’aire 


Approximation d’une aire par la méthode de Monte Carlo 


On souhaite de nouveau calculer l’aire sous la courbe de la fonction inverse entre 1 et 2. 
Ce domaine est contenu dans un carré ABDC que l’on a représenté ci-dessous. 
Figure. 


-0.24 


On choisit 10 points au hasard dans le carré ABDC : parmi ces dix points 6 se trouvent en-dessous 
de la courbe, alors que 4 points se trouvent au-dessus. 

La probabilité de choisir un point du carré sous la courbe est approximativement égale à 0,6 . 

Or cette probabilité est le rapport entre l’aire sous la courbe et l’aire du carré : une estimation de 
l’aire sous la courbe est donc 0,6 . 


Questions 


1. On veut améliorer la précision du résultat. 
On note n le nombre de points choisis au hasard dans le carré ABDC . 


1 
Un point M de coordonnées (x, y) est en-dessous de la courbe lorsque y < — 


Dans un premier temps, on considère l’algorithme suivant où T désigne le nombre de points en- 
dessous de la courbe : 
i ,n , X , y , I sont des nombres 
n prend la valeur 100 
POUR i allant de 1 à n 
x prend la valeur 1 + random () 
y prend la valeur random () 
Sy d RrALORS sestiiieureses,. 
FIN de POUR 
afficher ............... 
Compléter cet algorithme de telle façon qu’il affiche en sortie une estimation de l’aire sous la courbe. 
2. On souhaite maintenant que l’utilisateur indique en entrée le nombre n de points à choisir. 
Proposer et tester un algorithme qui réponde à cette question. 
Pour quelle valeur de n obtient-on expérimentalement une erreur inférieure à 1071 ? 
3. Déterminer une valeur approchée de x par la méthode de Monte Carlo. 


Corrigé Calcul d’aire (Monte Carlo) 


1. Première estimation de l’aire . 
i ,n , X , ÿy , I sont des nombres 
n prend la valeur 100 
POUR i allant de 1 à n 
x prend la valeur 1 + random () 
y prend la valeur random () 
SI y <= 1 / x ALORS T prend la valeur T + 1 
FIN de POUR 
afficher T / 100 
2. Estimation avec choix du nombre de points . 
i ,n , X , ÿy , I sont des nombres 


lire n 
EbC rise 
Résultats numériques. 
n 10 100 1000 10 000 
Pr 0,8 0,68 0,692 0,6983 


On peut estimer que si n > 100, l’erreur faite en approchant l’aire par cette méthode est inférieure 
à 1071. 
3. Il nous faut un domaine dont l’aire soit en rapport avec le nombre 7 : 
- le disque de centre O et de rayon 1 a pour aire x et pour équation x? + y? < 1 
- si on impose y > 0, on obtient un demi-disque de centre O et de rayon 1 qui a pour aire T/2 
- si on impose x > O0 et y > 0, on obtient un quart-de-disque d’aire 7/4 


La méthode de Monte Carlo est facile à adapter au cas du quart-de disque : on choisit au hasard des 
points dans le carré 0 < x < 1 et 0 < y < 1 et on teste pour chaque point si à? +y? <1. 
i ,n , X , y , I sont des nombres 
lire n 
POUR i allant de 1 à n 
x prend la valeur random () 
y prend la valeur random () 
SI xxx + y*y <= 1 ALORS T prend la valeur T + 1 
FIN de POUR 
afficher 4 x T/n 
Résultats numériques. 


n 10 100 1000 10 000 
A*T/n 2.4 3.08 3 172 3127 


Algorithmique n° 8 Simulation d’une expérience aléatoire 


Problème du Duc de Toscane 


On lance trois dés bien équilibrés. 
Est-il plus probable que la somme des chiffres obtenus soit égale à 9 ou à 10 ? 


Questions 


1. Montrer qu’il y a autant de façons de décomposer 9 ou 10 en une somme de trois entiers compris 
entre 1 et 6 , ces entiers étant pris dans un ordre quelconque. 
2. L’algorithme suivant simule l’expérience consistant à lancer trois dés et à calculer leur somme. 
i,n,X,Yÿ,Z,S , f et g sont des nombres 
n prend la valeur 1000 
POUR i allant de 1 à n 


x prend la valeur 1 + floor(6 * random ()) 
y prend la valeur 1 + floor(6 * random ()) 
z prend la valeur 1 + floor(6 * random ()) 
s prend la valeur x + y +2z 
SE ALORS. 5.5 ste 

FIN de POUR 

f prend la valeur g /n 

aTFiChÈT siemens 


Compléter l’algorithme de façon qu’on obtienne une approximation de la probabilité de l’événement 
« la somme est égale à 9 » 

3. Comparer le résultat précédent avec la probabilité de l'événement « la somme est égale à 10 » 

4. Est-il plus probable d'obtenir une somme égale à 9 ou à 10 ? Justifier votre réponse. 


Corrigé Simulation d’une expérience aléatoire 


1. Décompositions de 9 en une somme de trois entiers compris entre 1 et 6 : 
1 


DDNOR KH 
+ + + + + + 


3 
Décompositions de 10 en une somme de trois entiers compris entre 1 et 6 : 
1 


CO 
+ 
(en) 


(en) 
+ 
O 


1 

2 

2 

2 
3+3+4 

2. Approximation de la probabilité de l’événement « la somme est égale à 9 ». 


i,n,X,Yÿ,Z,S , f et g sont des nombres 
n prend la valeur 1000 


POUR i allant 
x prend la 
y prend la 
z prend la 
s prend la 


SI s = 9 ALORS g prend la valeur g + 1 


de 1 à 
valeur 
valeur 
valeur 
valeur 


n 
1 
1 
1 


X 


+ + + 


+ 


floor(6 * random ()) 
floor(6 * random ()) 
floor(6 * random ()) 


YF 


FIN de POUR 
f prend la valeur g /n 
afficher f 
Résultats numériques. 
n 10 100 1000 10 000 
g/n 0,1 0,2 0,122 0,1176 
. Probabilité de l’événement « la somme est égale à 10 » 
n 10 100 1000 10 000 
g/n 0,2 0,08 0,134 0,118 


4. Est-il plus probable d'obtenir une somme égale à 9 ou à 10 ? Difficile de conclure avec les données 
expérimentales. 
Il y a 6° = 216 issues toutes équiprobables. 
Une décomposition en trois entiers distincts est réalisée par 6 triplets. 
Si une décomposition comprend deux fois le même nombre, elle n’est plus réalisée que par 3 triplets. 
Enfin, si une décomposition comprend trois fois le même nombre, elle est réalisée par un seul triplet. 
La somme 9 est donc réalisée par 3 x 6 +2 x 3 +1 x 1 = 25 triplets . 
Alors que la somme 10 est réalisée par 3 x 6 +3 x 3 — 27 triplets. 
Il est donc plus probable que la somme soit 10. 


Algorithmique n° 9 Évolution d’un processus 


Le marché du téléphone 


Dans une ville imaginaire, le marché des télécommunications est partagé entre deux compagnies A et B. 
Initialement, la compagnie À truste 80% des abonnements. Mais la compagnie B est très agressive et 
cherche à augmenter sa part de marché. 
Dans cette ville, la règle veut que chaque habitant doit à la fin de chaque année confirmer son abonnement 
auprès de la compagnie qu’elle a choisie ou au contraire changer de compagnie. 
Une étude statistique conduit aux résultats suivants : 
- il y à 6 chances sur 10 qu’à la fin d’une année une personne cliente de À souscrive un abonnement 
auprès de B 
- il y a 3 chances sur 10 qu’à la fin d’une année une personne cliente de B souscrive un abonnement 
auprès de À 


Quelle compagnie aura à long terme la plus grosse part de marché dans cette ville ? 


Questions 


On note a, la proportion d'habitants qui sont clients de la compagnie À après n années. 
De même b,, désigne la proportion d'habitants qui sont clients de la compagnie B après n années. 
On a donc au départ ao = 0,8 et bo = 0, 2. 


1. L’algorithme suivant détermine l’évolution de la part de marché de chaque compagnie. 
i , n , àa , b , À , B sont des nombres 
a prend la valeur ......... 
b prend la valeur ......... 
n prend la valeur ......... 
POUR i allant de 1 à n 
À prend la valeur ........................ 
B prend la valeur ........................ 
afficher n 
afficher À 
afficher B 
a prend la valeur A 
b prend la valeur B 
FIN de la boucle POUR 
. Compléter et tester l’algorithme. Quelle conjecture peut être proposée ? 
3. a) Démontrer que ay+1 = 0, 4a» + 0, 3by 
b) En déduire ay+1 = 0, lan + 0,3 


D 


1 
c) Démontrer que un = an — 3 est géométrique 


d) En déduire la limite de (u,) puis de (a). 
4. Quelle réponse peut-on apporter à la question initiale ? 


Corrigé Evolution d’un processus 


1. D'après le texte : 
Ga 0,4an +0,30n 
bn+1 0,6 an + 0,7bn 
Évolution de la part de marché de chaque compagnie. 
i , n , àa , b , À , B sont des nombres 
a prend la valeur 0,8 
b prend la valeur 0,2 
n prend la valeur 10 
POUR i allant de 1 à n 
À prend la valeur 0.4 *x a + 0.3 x b 
B prend la valeur 0.6 *x a + 0.7 *x b 
afficher n 
afficher À 
afficher B 
a prend la valeur A 
b prend la valeur B 
FIN de la boucle POUR 
2. Résultats numériques. 


n 0 I 2 3 { D 
An 0,8 0,38 0,338 0,334 0,333 0,333 
bn 0,2 0,62 0,662 0,666 0,667 0,667 


1 2 
Conjecture : à terme les compagnies À et B auront respectivement 3 et — du marché. 


3. a) Les clients de À pendant l’année n + 1 sont constitués de 40% des clients de À pendant l’année 
n et de 30% des clients de B pendant l’année n. 


b) On remplace by par 1 — an : an+1 = 0,4an +0,3(1 — an) = 0,1 an + 0,3 


c) La suite auxiliaire un = an — 3 est géométrique : 


1 
Un+1 — An+1 — 3 


an 
10 I 


0, lun 


8 1 +4 
Sachant u9 = — — - — TE 


10 3 
7 LP 
SL I0 TS 


1 SALONS 
d) Puisque 1 < y < 1, on a lin (5) =; 


, on peut écrire Un = 


Gi 


On en déduit : an — 


On en déduit lim uy = 0 puis lim ay = a: 


4. La conjecture précédente est vérifiée. 


Algorithmique n° 10 Loi normale 


Approximation d’une loi binomiale par la loi normale 


Dans cet exercice, n désigne un entier naturel strictement positif. 
Au départ, on suppose n — 10 et on étudie l’expérience aléatoire suivante : 
- une pièce de monnaie est lancée n fois de suite 
- on compte le nombre de fois où on a obtenu « pile » au cours des n lancers 
Cette expérience aléatoire peut-être simulée par l’algorithme suivant : 
j ,; P et X sont des nombres 
X prend la valeur 0 
POUR j allant de 1 à 10 
P prend la valeur floor(2 *xrandom()) 
SI P = © ALORS X prend la valeur X + 1 
FIN de la boucle POUR 
Afficher X 
Dans cet algorithme, le lancer de la pièce est simulé par le choix d’un nombre entier aléatoire égal à 0 
(pile) ou à 1 (face). 
On note désormais X la variable aléatoire égale au nombre de « pile » au cours d’une série de n lancers. 


Questions 


1. Quelle loi de probabilité est suivie par X ? 
En déduire l’espérance y, la variance V et l’écart-type o de X. 
2. On veut déterminer expérimentalement p(X = 5) en simulant un grand nombre de fois l'expérience 
initiale. Proposer et tester un algorithme qui réponde à la question. 
On pourra utiliser les variables suivantes : 
N : nombre de simulations (déterminé en entrée par l'utilisateur) 
i : compteur associé aux simulations (ce compteur est distinct du compteur j) 
G : nombre de cas où X prend la valeur 5 
F : fréquence associée à l'événement (X = 5) qui doit être affichée en sortie 
Indiquer les valeurs obtenues dans un tableau et comparer avec la valeur théorique. 
3. Montrer que l'événement (u — o < X < u +0) est la réunion des événements (X = 4), (X —5) et 
(CE = 6): 
Adapter l’algorithme précédent pour obtenir une approximation de p, = plu —-o<X<p+o). 
Tester et comparer avec la valeur théorique. 
4. Dans cette question, on va estimer p1 par une méthode différente. 


Sur cette figure chaque rectangle a une aire égale à p(X = k) où k est l’abscisse indiquée en-dessous 
du rectangle. 


On a également représenté une fonction f qui « approche » le dessin constitué par ces rectangles : 
cette fonction s’appelle « la densité de probabilité de la loi normale d'espérance 1 = 5 et de variance 
o? = 2,9 ». 

On démontre qu’une très bonne approximation de p1 est l’aire sous la courbe de f entre les abscisses 
a=yp—oetb=u+a.On va déterminer l’aire de ce domaine avec différents outils. 

Dans chacun des cas considérés, les variables a , b , u et « devront être remplacées par leurs valeurs 
numériques. 


a) Lancer Geogebra et définir dans la barre de saisie la fonction densité : 
f(x) = Normalelu , o , x] 
Déterminer l’aire sous la courbe entre les bornes indiquées : 
IntégraleDomainel[f , O , a , b] 


b) Lancer Algobox et définir la fonction numérique : 
F1(x) = ALGOBOX LOI NORMALE(U , © , x) 
Cette fonction calcule l’aire sous la courbe de f qui est située à gauche de l’abscisse x. 
Il suffit donc dans un algorithme de calculer la probabilité recherchée en écrivant l’instruction : 
pi prend la valeur F1(b) - F1(a) 


c) Lancer un tableur et utiliser la fonction LOI . NORMALE 
Cette fonction renvoie la valeur de la densité en x avec : 
= LOI.NORMALE(x ; H ; 0) 
La même fonction renvoie l’aire sous la courbe de f qui est située à gauche de l’abscisse x avec : 
= LOI.NORMALE(x ; j ; © ; 1) 
Il suffit donc de taper dans une cellule du tableur : 
= LOI.NORMALE(Db ; y ; © ; 1) - LOI.NORMALE(a ; j ; © ; 1) 


d) Vérifier les résultats obtenus avec votre calculatrice 
TT : utiliser la fonction normalcdf(a , b , H , o) 
CASIO : dans le menu STAT DIST NORM Ncd préciser les valeurs de Lower , Upper , o et y 


. Application. On lance 100 fois une pièce de monnaie et on note X le nombre de « pile ». 
En convenant de noter u l’espérance de X et o son écart-type, déterminer la probabilité des événe- 
ments suivants : 

(X < 40) 

(X > 90) 
(H=TSX< io) 
(u — 20 < X < u+20) 
(u — 30 < X < u +30) 


Corrigé Approximation d’une loi binomiale 


1. X suit une loi binomiale de paramètres n = 10 et p = 0,5 
L’espérance est u = np = 5, la variance est V = np(1—p) = 2,5 et l’écart-type est o = V2 & 1,58 
2. Estimation de p(X = 5) 
i,j, P,F,G, Net X sont des nombres 
lire N 
POUR i allant de 1 à N 
X prend la valeur 0 
POUR j allant de 1 à 10 
P prend la valeur floor(2 *xrandom()) 
SI P = O ALORS X prend la valeur X + 1 
FIN de la boucle POUR 
SI X = 5 ALORS G prend la valeur G + 1 
FIN de la boucle POUR 
F prend la valeur G / N 


Afficher F 
Résultats numériques pour l'événement (X = 5). 
N 10 100 1 000 10 000 
F 0,3 0,231 0,239 0,2495 


10 
Remarque : la valeur théorique est p(X = 5) — ( 5 L 510 & 0,246 


3. On vérifie p — o & 3,42 et p + o & 6,58. Par conséquent : 
p—o<X<up+o — X =4 où X =5 ou X —=6 
Algorithme permettant d'estimer p; = p(E —-o < X < E +0) : il suffit de remplacer dans l’algo- 
rithme précédent le test 
X = 5 
par le test : 
4 <= X ET X <= 6 

Résultats numériques pour p1 : 

N 10 100 1 000 10 000 


F 0,5 0,61 0,668 0,6574 


Calcul de la valeur théorique : 


pi = p(X = 4) + p(X = 5) + p(X = 6) 


— 0,5 0,5 0,5 
() +() 5+(, 0 


= 672 .0,510 & 0,656 
4. Approximation par la loi normale. 
a) Geogebra 


b) Algobox 


p prend la valeur F1(5 + 1.58) - F1(5 - 1.58) 
afficher p 


Fonction numérique utilisée : F1(x) = ALGOBOX LOI NORMALE(5 , 1.58 , x) 


c) Tableur : OpenOffice 
= LOI.NORMALE(6,58 ; 5 ; 1,58 ; 1) - LOI.NORMALE(3,42 ; 5 ; 1,58 ; 1) 


d) Calculatrices 
TT: normalcdf(3.42 , 6.58 , 5 , 1.58) 
CASIO : renseigner Lower Upper o& y dans le menu Normal C.D. et calculer 


. Application au cas de 100 lancers. Dans ce cas : 
l'espérance est 1 = 100 X 0,5 = 50 
la variance est V — 100 x 0,5 x 0,5 = 25 
l’écart-type est o = V25 = 5 

Les résultats suivants ont été obtenus avec un tableur : 
p(X < 40) = 0,02275 


